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Stellar wobble caused by a binary system: 
Can it really be mistaken as an extra-solar planet? 
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ABSTRACT 



The traditional method for detecting extra-solar planets relies on measuring a small stellar wobble which is assumed to be caused by 
a planet orbiting the star. Recently, it has been suggested that a similar stellar wobble could be caused by a close binary system. Here 
we show that, although the effect of a close binary system can at first sight be mistaken as a planetary companion to the star, more 
careful analysis of the observational data should allow us to distinguish between the two effects. 
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ed 1. Introduction 



Over the past decade up to 290 extra-solar planets were dis- 
covered using the wobble method (i.e. measuring the motion of 
parent star around the planet-star's center of mass). This repre- 
sents about 95% of the total extra-solar planets known to dateQ 
However, this is an indirect method in which the presence of a 
planet is inferred from the detection of the radial velocity vari- 
ations. Therefore, in cases where the presence of the planet can 
not be confirmed by other methods we may ask if these vari- 
ations could be due to another dynamical effect. In particular, 
when a planet is detected in a system already hosting a compan- 
ion (such as another planet or a star), we must be cautious when 
we analyze the data, because different configurations of bodies 
can lead to similar radial velocity variations. 

The present pa per is motivated by the work of 

■ ISchneider & Cabreral ([2006) who studied a triple system 
composed of a binary star system with equal masses orbiting a 
third star (assumed massless). They concluded that the binary 
system will cause a wobble in the star's motion that could 
mimic the presence of a planet companion. Therefore, they 
suggest that each extra-solar planet detection which relies on 
the measurement of this wobble must be carefully checked for 

the presence of possible nearby u nresolved binaries. 

Unfortunately the study in ISchneider & Cabreral d2006l) was 
restricted to the case of binary systems composed of stars with 
equal masses, while a multitude of three-body configurations can 
be imagined. In particular, we may wonder about the effect of 
a giant satellite orbiting an already detect ed Jupiter-like extra- 
solar p lanet. Moreover, the conclusions in ISchneider & Ca brera 
(2006) were derived in the framework of the restricted three- 
body problem and we may also question the validity of this 
approximation. Additionally, we found that their work contains 
several inaccurate conclusions. 
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Our goal is to correctly model the effect of a binary system 
on a nearby star without making any assumptions on the masses 
of this triple system. In Sect. 2 we briefly review the wobble 
method for detecting planets. In Sect. 3 we derive simple ex- 
pressions for the radial velocity of a star under the presence of a 
nearby binary system and for the orbital parameters of a planet 
mimicked by this effect. In Sect. 4 we give some examples of bi- 
nary systems and present numerical simulations to test the theo- 
retical results. Finally, the last section is devoted to a discussion 
of the results. 

2. The wobble method 

We quickly review the wobble method for detecting extra-solar 
planets (see for instance iHildi tch 2001). Assume that we are ob- 
serving a star with mass M*, which has a planet companion with 
mass M p . A possible observation reference frame is depicted in 
FigQJ where we have the x and z axes along the line of nodes 
and the line of sight, respectively. The orbital plane (i, j) is de- 
picted in FigOJ), where i = x. In this last frame the position of 
the star/planet can be written as 



r — X i + Yj — r cos(-nr + /) i + r im(vj + f)j 



(1) 



where r is the orbital radius, / the true anomaly, and w the lon- 
gitude of the perihelium. 

Since i = x, the j axis is simply obtained with a rotation 
of the y axis by an angle /, the inclination of the orbital plane 
measured with respect to the plane (x,y). Thus, we may write 
z = Y sin /, where z is the position along the z axis. The velocity 
along the line of sight is then 



Fsin/ 

(r sin(cr + f) + rf cos(?zr + /)) sin / 
K z [cos(nr + f) + e cos m] , 



(2) 



where e is the orbital eccentricity and K z = na p sinl/ Vl - e 2 , 
with n the mean motion, and a p the semi-major axis of the orbit. 
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Fig. 1. Reference frames. The plane defined by (x,y) is the plane 
of the sky and the z axis is along the line of sight. The plane (i, j) 
is the orbital plane, m the argument of the perihelium and / the 
true anomaly. For simplicity we chose the x axis along the line 
of nodes (i = x). 



The radial velocity of the star is the projection of the velocity 
with respect to the center of mass (CM) along the line of sight. In 
the particular case of a circular orbit, we have e — and f - n, 
thus 



V r 



Mn 



■ z = K cos( m + nt) . 



where t is the time and the amplitude 



K 



Mn 



M++M, 



n a p sin / . 



(3) 



(4) 



If we analyze the Fourier spectrum of the star's radial ve- 
locity we will see that this will exhibit a peak at the frequency n 
and the amplitude K provides a lower limit for the planet's mass. 
This method of detection of extra-solar planets is also known as 
the wobble method (since the planet causes the star to orbit or 
wobble about their common CM). 

Notice also that when more than one companion to the star 
is present (another planet or star), the total radial velocity will 
be given by a sum of the individual contribution from each body, 
V r = ZkM k Zk/M Tot , where M Tor = M+ + ZkM k is the total 
mass of the system. The amplitude due to a specific companion 
is then: 



M p 

K n — — na„ sin / 

' M*+X k M k p 



(5) 



3. The model and its predictions 

We will work in the framework of the planar general three body 
problem: we assume that a star with mass M* moves in the same 
orbital plane as the binary system composed of two bodies with 
masses M\ and M 2 (Fig. [2j. 




Fig. 2. Jacobi coordinates (r b ,r). The system is composed by a 
star with mass M* and a binary with masses M\ and Mi. We 
assume co-planar motion and the frame (i, j) is the same from 
FigQ] u is a unit vector normal to r, S = L(r, r b ) and 9 = r). 



3.1. Equations of motion 

In order to obtain the equations of motion we will use Jacobi 
canonical coordinates, which are, respectively, the inter-binary 
distance, r B , and the distance, r, from the star M* to the 
binary's CM (see Fig. |2). We will additionaly assume that 
\r B \ «: \r\. Thus, to the second order in p — \r B \/\r\ < g 1, the 
Hamiltonian of the system can be expressed as (e.g. iHenrardl 
ll979HMurray & Dermottlll999h : 



H = 



1W 2 1\P\ 2 n M,M 2 MM* 

H O O 

2 ju* r b r 



2 

-G 



M 1 M 2 M i , 



(6) 



p 2 l(3cos 2 S-l) 



rM r 2 

where M - M\+ M 2 is the total mass of the binary, S - l(r, r D ), 
M\M 2 , MM+ 



V = 



M 



and 



(7) 



M + M+ 

are the reduced masses and p D = firb and p = n+r the conjugate 
momenta of r and r, respectively. In this framework, since the 
Jacobi coordinates are canonical, the equations of motion are 
easily obtained as 



dH , dH 

Pb = -3— and P = — r- • 
or b or 



3.2. Motion of the binary 

From Eq.© the equation of motion for r b is 

M M* 
ir b = -G-rr b + G— -(3p cos Sr - r b ) 

rt 



(8) 



(9) 



We can neglect the second term in Eq.© if p 3 <K M/M+, which 
is always true if we choose a "stable binary", i.e. if Mj > M 2 , 
such that 



Mi 

3m7 



1/3 



(10) 
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where a 5 0.5 ( iMarkellos & Rovll981l) . Therefore, the 0th order 
solution for rt is a Keplerian ellipse with constant semi-major 
axis at, and frequency 



n/co 



jt/Q 



LO - 




(ID 



3.3. Motion of the star 

From Eq.© the equation of motion for r is 



(M + M*) 
(l + ^y(-3+15cos 2 5)]r 



(12) 



M 



3p cos Srt 



We can neglect the 2nd and 3rd terms in Eq.dT2b if 
« 1 . 



(13) 



Since p/M < 1/4 and from Eq.dTOt p 2 <K 1, this is always true. 
Therefore, the 0th order solution for r is also a Keplerian ellipse 
with constant semi-major axis a and frequency 



Q. 



G(M + M*) 



(14) 



i.e., the motion of M* is to a first approximation a Keplerian 
orbit about a mass M located at the binary's CM. Furthermore, 
from Eqs. ( fTTT i and < TT~4T > we see that in general Q. <*c a> since 
at, <K a (this is true unless M* becomes too large). 

3.4. Stellar wobble 

Now, we can obtain an approximation for the relative motion of 
M* by replacing the 0th order approximations for rj and r in 
Eq.([T2l. We assume that these 0th order solutions are circular 
orbits, thus 



r — ar , 

rt = a/, cos Sr + at, sin 5 m 
5 = 5 + (w-Q)f , 



(15) 
(16) 
(17) 



where f is the versor of r, u is the unit vector orthogonal to r and 
5 o is an initial phase (Figj2j. 

In order to compute the radial velocity, we are interested in 
obtaining the equation of motion for M* in barycentric coordi- 
nates. The distance of the star to the CM is 



M 



r* = 



M + M+ 
We hence have 



SI- 



1 + — p 2 (- + - cos 25)1 f 
AT 4 4 1 



u j 5 . ^ _ „ 
— p — sin 25 m 



(18) 



(19) 



with p = a/, /a <k 1. Notice that if we chose M* = and Mi = 
Mo (binary with equal masses) in Eq.(fT9|> we recover the same 
equation of motion obtained by Schneider & Cabrera (2006). 




3co6 



50000 



52000 



54000 56000 
JD-2400000 (day) 



58000 



60000 



Fig. 3. Radial velocity of a star due to the binary wobble around 
its own CM (Eql25ll. The orbital parameters are from Ex.#3, 
Tabg] 



In Appendix A we determine the solution to Eq.([T9l in the 
inertial frame (i, j) with coordinates = (X, Y), that is: 



a*cos0\, 
a* sin 6 J ' 



sin 6 cos 6 
- cos 6 sin 



(3 5 sin 25 
4.5 8 cos 25 



with a* * aM/(M + M+), = O + Q, t and 

\4 



8M \ a ) 



a b . 



(20) 



(21) 



This solution has two components, the first one caused by the 
motion of the star around the binary's CM and the second one is 
the stellar wobble caused by the the movement of the two bodies 
of the binary around its own CM. 

According to our choice of reference frame (FigQ}, the radial 
velocity of the star M* is given by V, = Y sin /, that is, we can 
obtain it directly from the Y coordinate. From Eq.(l20l we have: 



Y — a* sin 6 - 3 5 sin 25 cos 6 + 4.5 5 cos 25 sin , 



(22) 



Because 5 = u> - Q. and 
the velocity along y is then 



Q, and since we assume £2 <K a>, 



Y = a+Q. cos - 6a>5 cos 25 cos - 9a>8 sin 25 sin . 



(23) 



The first term in Eq.(|23b is due to the slow motion of the star 
around the binary's CM. The second and third terms are due to 
the binary wobble. These last two terms can be seen as a compo- 
sition of two periodic signals, one fast with frequency 2co - 2Q. 
and another much slower with frequency Q. The global effect of 
the binary wobble corresponds to a signal of period 7i/((o - Q.) 
with an amplitude modulation of 3a>6 and period tt/O (FigfJJ). 

3.5. Mimicking a planet 

Even though the radial velocity of a star perturbed by a binary 
system is not given by a single periodic signal (Fig[3]), we may 
ask under what circumstances this signal can be misinterpreted 
as a planet orbiting the star. We can rewrite Eq.(f23]) in the same 
format as Eq.(|3j, i.e. 



V r = K o cos(0o + Q.t) + V° , 

where Kq - a+Q. sin / and 

V* = K\ cos(tzn + n\t) + K2 cos(tzr 2 + H2O , 



(24) 



(25) 
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Table 1. Planets detected in close binary systems and the respective parameters (at,, Mq) of the stellar companion that can mimic 
them. 



system name 


M, (M Q ) 


M, (Mo) 


M p (M Jup ) 


fl(AU) 


a P (AU) 


M 2 (M Q ) 


a b (AU) 


a 


GJ86A 1 


0.79 


-0.5 


4.01 


-18 


0.11 


3.4 x 10 3 


2.85 


0.27 


yCepA 2 


1.40 


0.41 


1.60 


20.2 


2.04 


3.81 


4.68 


0.50 


HD 41004 A 3 


0.7 


0.4 


2.54 


23 


1.31 


11.9 


5.41 


0.41 


HD 41004 B 3 


0.4 


0.7 


18.37 


23 


0.018 


4.9 x 10 5 


3.06 


0.16 


HD 196885 A 4 


1.33 


-0.6 


2.96 


-23 


2.63 


2.93 


5.78 


0.47 


HD 19994 A 5 


1.34 


0.3 


1.68 


-100 


1.42 


6.2 x 10 2 


17.4 


0.41 



Reference s: llllOueloz et alJ(|2000h:lLagrange et all j2006ft ; mlNeuhauser et alj d2007l) ; 
[3] IZucker et alj J2004h : MCorreia et alj d2008l) . r51 iMavor et al.1 d2004l) . 



with 

15 3 
K\ = -—codsml and Kz = —u>5sml. 

n\ - 2oj - 3Q and «2 = 2co - Q , 
and 

vj i = 25 o - 0q and mq — 2S o + 6>o . 



(26) 
(27) 

(28) 



The 1st term in Eq.(l24li represents the slow two-body motion 
around a body with mass M located at the binary's CM. We re- 
call that we may not know about the presence of the binary (since 
one or even both components may be unresolved) and this term 
allows us to identify only a single component. 

On the other hand, V* is the radial velocity due to the bi- 
nary wobble around its own CM and can be identified as two 
planets on circular orbits ar ound the star. Hence, contrarily to 
ISchneider & Cabreral ([2006), we conclude that the stellar wob- 
ble caused by the presence of the binary mimics not one, but 
two planets, although with very close orbital periods. As the two 
frequencies n\ and «2 are very close, in some situations the two 
signals can be confounded and the observer can erroneously be- 
lieve to have found a planetary companion to the star. 

This will be the case if our instrument is at the limit of its 
resolution: since \K\\ -5 \K 2 \, the trace of the second "planet" 
is harder to detect. Therefore, if our spectral resolution is close 
to the detection limit of the stellar wobble, we can notice the 
presence of the first "planet", but we will miss the second one. 

By applying Kepler's 3rd law, we obtain for the semi-major 
axis of this fake planet, 



(GM+) 



1/3 



(-) 

Kami 



1/3 



a h 



(29) 



and for the planet's mass, we replace K\ from Eq.d26li in Eq.© 
M„ 15 6 co 15 n 



M 



M p _ 15 b to ^ 15 n la h ^ a b 
+ M+ 2 a p n\ 32 M \ a I a p 



(30) 



Notice that when Mi = M 2 our estimate f or the planet's mass 
is diff erent from what was obtained by I Schneider & Cabreral 
(2006). Additionally, we also conclude that the planet evolves 
in a circular orbit (wh ile e = 0.745 from Eq.4 of 
ISchneider & Cabrera 20061) . 



3.6. Distinguishing between a binary and a planet 

The situation described above is very likely and we may ex- 
pect to interpret erroneously the signal of a binary as a planet. 
However, the radial velocity wobble produced by a binary 



(Eqf25]), is different from the radial velocity produced by a sin- 
gle planet (Eq{3j and therefore we may be able to distinguish 
between the two situations. 

Indeed, if we acquire data for long enough time (so that we 
can constrain well Q.) and if the precision of the instrument is 
at the limit of detection of the amplitude K 2 , the observer will 
be able to notice the presence of the two periodic signals «i and 
«2 (Eql27l> and therefore realize that the star is undergoing the 
perturbation from a binary and not from a planet. 

Since the radial velocity of the binary (Eqj25t is the same 
as the signal produced by two planets in circular orbits, the ob- 
server can also believe that two planets have been found (instead 
of the binary). However, as long as fl « u, the two orbital fre- 
quencies ni and «2 are very close and probably the orbits of the 
two planets cannot be stable. 

For higher £2 values, the system may become stable (al- 
though the approximation !) « w may not be valid) and 
the hypothesis of existence of two planets cannot be excluded. 
Nevertheless, in the case of a binary, the frequencies (m, «2) and 
the amplitudes (K\, K2) are not independent (Eqsf27]and|26l) and 
the following relations must be verified: 



n 2 -«i=2Q and \Ki\=5\K 2 \ 



3. 7. What kind of binary can mimic a planet? 



(31) 



In Sect. 13. 51 we saw that under some particular conditions, the 
variations in the star's radial velocity resulting from the binary 
wobble can be misinterpreted as a planet. Thus, if we detect a 
planet and we do not know about the presence of the binary sys- 
tem (e.g. if M2 is unresolved) we can ask the following ques- 
tion: what are the parameters of an hypothetical binary system 
that can mimic this planet? More precisely, if we know M+, Mi, 
a, M p and a p , what are the values of M2 and a^,? These can be 
obtained by inverting Eqs.(|29l and ( l30l >: 



MS 



(GM) 



1/3 



4M 

m7 



■('♦£) 



M„ 



(In, + §£2) 2/3 

16 n\ a p I a\ 
15 lo at, \at, J p 



1/3 



(32) 
(33) 



M„ 



In Table[T]we list all the currently known planets in close 
binary systems (a < 100 AU) and estimate the corresponding 
companions' parameters (M2 and at,) that can mimic that planet. 
Three of the cases provide unrealistic binary systems with very 
high values of M2. The remaining three casesare are acceptable 
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(y Cep A, HD 41004 A, and HD 196885 A) but these binary sys- 
tems are likely to be unstable as a is near the stability limit (see 
Eg HOIi. On the other hand, the values of M 2 in these three cases 
are in the mass range of either black-holes or very bright stars 
(which should be detectable). Therefore, it is unlikely that these 
planets can be the result of a binary wobble. 

4. Examples 

In Sect.|3]we obtained the radial velocity wobble of a star that is 
perturbed by a binary system. We subsequently derived theoreti- 
cal expressions that allow us to determine the orbital parameters 
of a planet that can, in some circumstances, be mimicked by 
those variations. In this section we will present some concrete 
examples of binary systems together with the orbital parameters 
of the planets that they can mimic. We also obtain the amplitude 
of the peaks in the radial velocity curve since these give an in- 
dication of the instrument precision needed to detect the planet 
and to distinguish between the planet and binary. 

4.1. Theory versus numerical integrations 

In order to test the theoretical predictions of our model (Sect.[3]l, 
we performed some numerical simulations of three-body sys- 
tems (Tab|4]i. Here we present in detail a case (Ex. #3, Tab|4|i for 
which the current precision in the radial velocity measurements 
(~ 1 m/s) could lead to the erroneous announcement of a planet. 

The hypothetical system is formed by a triplet of main se- 
quence stars with different spectral types, G, K and M, and 
masses M* = 1.00 M , Mi = 0.70 M and M 2 = 0.35 M@, 
respectively. The smaller K and M stars form a binary system on 
a circular orbit with flj, = 1.1 AU. The G-star is in a wider cir- 
cular orbit around the close binary's CM, with a = 10 AU. Note 
that the M star is much fainter than the K or G stars hence it 
represents the unresolved component of the binary. We also sup- 
pose that the system is co-planar and perpendicular to the line of 
sight (that is, sin / = 1). 

We numerically integrated this stellar system and computed 
the radial velocity of the G-star, the brightest body in the system. 
In Table [2] we list 50 simulated observational data points for a 
time span of about 10.5 years (from Feb. 2004 to Sep. 2014). We 
supposed that the data was acquired with a precision of ~0.8 m/s, 
identical to the HARPS spectrograph on the ESO 3.6-m telescope 
at La Silla (Chile), currently the mo st precise instrum ent in op- 
eration since the beginning of 2004 dMavor et al.ll2003l) . 

From Eq.dTTI) and Eq.(fT4b we compute for this system: 

— = 411.23 day and — = 8067.0 day . (34) 

ID U. 

The observational data do not fully cover the orbital period of the 
G-star around the binary's CM, but we will be nevertheless able 
to estimate it, since we are assuming a high level of precision in 
the data. 

According to our model predictions (Sect. l3.5l l. the radial ve- 
locity of the G-star can be decomposed in three sinusoidal terms 
(Eql24li. One of large amplitude, Kq = 6.907 km/s, but varying 
slowly with period 2n/Q, and the other two terms varying much 
faster with periods (EqlZTjl 

2n 2n 



— = 222.64 day and — = 

«1 ' «2 

and small amplitudes (Eql26t 



210.99 day, (35) 



|#i| = 0.89 m/s and \K 2 \ = 0.18 m/s . (36) 



Table 2. Simulated radial-velocity data for a triple system of 
stars with M* = 1.00 M , M x = 0.70 M , M 2 = 0.35 M B , 
a = 10.0 AU and ah - 1.1 AU. The 50 points are generated for a 
period of about 10.5 yr. 
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Since the precision of the instrument is supposed to be around 
0.8 m/s, it should be able to detect the signal from Ky, but not 
from K 2 . 

With the 50 data points listed in Table |2] we will now apply 
the traditional techniques used to search for planets. The orbit 
of the binary CM introduces large amplitude variations in the 
radial velocity of the G-star and at first glance we can only de- 
tect thes e variations. Us ing the iterative Levenberg-Marquardt 
method (iPress et al.ll 1992b . we then first attempt to fit the com- 
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Fig. 4. Radial velocities for the G-star superimposed on a 
Keplerian orbital solution due to the close binary CM. The re- 
maining residuals are slightly above the precision of the instru- 
ment and the observer can suspect of the existence of second 
companion in the system. 




1 10 100 1000 10000 

Period [days] 



Fig. 5. Frequency analysis (top) and periodogram of the dates 
(bottom) for residual radial velocities when the contributions 
from the binary CM is subtracted. An important peak is detected 
at P — 223 day, which can be interpreted as a second body in the 
system. 



plete set of radial velocities with a single orbiting companion in 
a Keplerian orbit. This fit yields a stellar companion with P = 
8010 day, e = 0.002 and a minimum mass of 1.044 M Q (FigJUi. 
These parameters correspond very well to the orbit of the G-star 
around the binary's CM. The mass determined by the fit is the 
total mass of the close binary system, M — M\+ Mi = 1.05 M Q . 

In Figure H] we plot the radial velocity of the star superim- 
posed on the binary CM orbit and also the velocity residuals after 
subtracting these long term variations. The precision of the in- 
strument is ~1 m/s, but the residuals, although very small, show 
amplitude variations that can reach almost +3 m/s. A trained ob- 
server can then suspect for the existence of something else in the 
data. To check it in a simple way we perform a frequency anal- 
ysis of the residuals (Fig|5). We can then see an important peak 
around 223 day, clearly indicating that something else is hidden 
in the data. We also plot a periodogram of the data, and since 
there is no peak at the same frequency (there is only a small one 
around 364 day) we conclude that the signal is true. 



Table 3. Orbital parameters for two bodies orbiting the G-star, 
obtained with a two-Keplerian fit to the observational data. 



Param. 


[unit] 


planet 


binary CM 


i/ 
Vq 


LKin/sj 


i n tv> 1 

1U.UZ1 


± U.U9 / 


P 


|uayj 


999 f.] + 1 99 


OUIO.U It Ji.J 


e 




0.0486 ±0.1725 


0.0016 ±0.0002 


CO 


[deg] 


94.04 ± 237.2 


280.35 ± 34.93 


K 


[m/s] 


1.07 ±0.19 


6904.6 ±30.1 


T 


[JD-2400000] 


53839.0 ± 146.9 


53471.9 ± 778.0 


a\ sin i 


[10- 3 AU] 


0.022 


5089 


f(m) 


[io- 3 M Q ] 


2.8 x 10-" 


273.5 


m.2 sin i 


[M e ] 


20.6 


3.5 x 10 5 


mt sin f 


[Mo] 


6.2 x 10~ 5 


1.046 


a 


[AU] 


0.719 


9.95 


rms 


[m/s] 


0.56 


V? 




0.51 



Errors are given by the standard deviation cr. The orbital period of the 
outer body is longer than the data acquired and thus we are unable to 
complete constrain its orbit. The eccentricity and the perihelium of the 
inner planet are also very badly constrained because we are at the limit 
of the instrument resolution. Despite all the uncertainties, the adjusted 
parameters for the planet agree with our model predictions CSect. l3.5l >. 



We then fit the radial velocities (Tabf2|i using a model with 
two Keplerian orbits. It yields for the planetary companion Pi = 
222.6 day, e = 0.049 and a minimum mass of 20.6 M m , while 
for the binary's CM we get approximately the same orbital pa- 
rameters as before (Tab|3). Despite all the uncertainties in the 
orbital parameters, the use of a second Keplerian orbit is a better 
approximation than the single Keplerian model, since the veloc- 
ity residuals drop to rms = 0.56 m/s, below the precision of the 
instrument (they were above 1 m/s before). 

In Figure|6]we plot the phase-folded residual radial velocities 
when the contribution of the binary CM is subtracted. The data 
is superimposed on the orbital solution of the planetary compan- 
ion, which really looks like a planet. The remaining O - C is be- 
low the precision of the instrument (~ 1 m/s) and therefore we do 
not find any additional signal in the data. This is confirmed when 
we perform a frequency analysis of these residuals: no relevant 
peak is found, and there is no signal around P2 = 211 day, the 
period for which we expected to see the «2 frequency (Eqj35j. 

By comparing the numerical determination period of the 
planet, Pi = 222.61 ± 1.22 day (Tab|3]l with its theoretical es- 
timation Pi = 222.64 day (Eq 13511. we see that our prediction 
was correct. Results for the amplitude are not so good, but still 
exact since the numerical determination K\ = 1.07 ± 0.19 m/s 
(Tabj3]l includes within the error bar the theoretical estimation 
K\ = 0.89 m/s (Eq|55l. A better numerical determination could 
be achieved with identical instrumental precision, but using a 
larger amount of data points. On the other hand, the amplitude 
K 2 = 0.18 m/s (Eqfjg) is lower than cr Kl = 0.19 m/s (TabfJ), the 
error bar of K\ , so there is no hope of determining it with the 
present experience. 

From the numerical values of n\ and K\ we can infer the 
orbital parameters of the planet, a p = 0.719 ± 0.003 AU and 
M p = 20.6 + 3.7 M e , while the theoretical values were respec- 
tively a p = 0.719 AU (Eqj29j» and M p = 17.3 M m (EqUJ. The 
mass is overestimated because the amplitude K\ was so. 

Now that we know that a binary can mimic a planet, when- 
ever we suspect that a binary can be responsible for the observed 
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Table 4. Model predictions for different examples of binary systems around a star with M+ = M Q . 



1 







Binary system 






Frequencies 


Amplitudes 


Planet 


Ex. 


My 


M 2 


a 


a b 


a 


2n/Q. 


2n/tii 




\K X \ 


\K 2 \ 


M p 


a p 




(M ) 


(M Q ) 


(AU) 


(AU) 




(yr) 


(day) 


(day) 


(m/s) 


(m/s) 


(M ffi ) 


(AU) 


1 


1.00 


1.00 


10.0 


1.00 


0.114 


18.26 


137.1 


131.7 


0.987 


0.197 


23.90 


0.520 


2 


1.00 


1.00 


10.0 


1.50 


0.172 


18.26 


265.6 


246.0 


4.081 


0.816 


123.1 


0.809 


3 


0.70 


0.35 


10.0 


1.10 


0.156 


22.09 


222.6 


211.0 


0.888 


0.178 


17.26 


0.719 


4 


1.00 


0.10 


10.0 


1.50 


0.210 


21.82 


363.7 


333.3 


1.000 


0.200 


23.47 


0.997 


5 


1.00 


0.01 


10.0 


1.50 


0.215 


22.30 


380.6 


348.1 


0.114 


0.023 


2.59 


1.028 


6 


1.00 


IO" 3 


10.0 


1.50 


0.216 


22.36 


382.5 


349.7 


0.012 


0.002 


0.26 


1.031 


7 


io- 3 


io- 3 


1.00 


0.01 


0.114 


1.00 


4.23 


4.13 


io- 5 


10~ 6 


io- 4 


0.051 


8 


io- 3 


10~ 6 


1.00 


0.01 


0.144 


1.00 


6.06 


5.86 


io- 7 


io- 8 


io- 7 


0.065 
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Fig. 6. Phase-folded residual radial velocities when the contribu- 
tion of the binary CM is subtracted. Data is superimposed on a 
Keplerian solution with P = 222.6 day and e » 0. The respective 
residuals as a function of Julian Date are displayed in the lower 
panel. We see that the remaining O - C is below the precision of 
the instrument (~1 m/s) and therefore we do not expect to find 
any additional signal in the data. 



stellar wobble, we should also provide its parameters. From the 
values of n\ and Q. (Tabj3) we can obtain (Eq|27t 



In 

0} 



41 1.0 ±2.2 day and 



2tt 



8018 + 33 day, 



(37) 



which are in very good agreement with the true values (Eql34l). 
The individual masses of the stars in the close binary can also be 
resolved as well as the semi-major axis. Directly from the fit we 
have a = 9.95 + 0.03 AU and M = 1.046 + 0.005 M Q (Tab©. 
Then, using Eqs.(l32l and d33l we can respectively obtain 



a b = 1 .098 + 0.006 AU and yu = 0.275 + 0.062 M e 



(38) 



By comparing the above determinations with the true pa- 
rameters chosen for this system (Ex. #3, Tab© we verify that 
they match quite well. We thus conclude that a two-Keplerian 
approach of the stellar radial velocity is enough to determine 
the parameters of the triple system, i.e., we can keep the current 
software that is used to fit planets to the stellar data and do not 
need to develop a specific algorithm for binaries. Of course this 
in only true if we keep a Keplerian approach, the use of a N-body 



algorithm will give different results in the case of a planet or a 
binary. 

4.2. Different kinds of binary systems 

Now that we have seen that our formulae are correct we will 
apply them to a wide variety of binary systems' configurations. 
In Table|4]we see the predictions of the model for several binary 
systems around a star with M* = Mq. We chose only stable 
binary systems with a « 0.5 (see EqflOb. 

The first and second e xamples are the cases studied in 
ISchneider & Cabreral ((2006): the star is being perturbed by close 
equal masses' binaries composed of stars with 1 Mq. In Ex.#l, 
K\ is at the current resolution level (~1 m/s) while is 5 times 
smaller. In this case we may be lead to believe that the star has a 
planet companion with 1.4 Neptune masses. On the other hand, 
in Ex.#2 we should be able to detect both K\ and K%\ hence in 
this case we can conclude that the star's wobble is caused by a 
binary and not a planet companion. 

The third example is the case discussed in the previous sec- 
tion. Here we have a triple system composed of stars with dif- 
ferent spectral types. As we saw previously, in this case we can 
detect the peak with amplitude K\ , but not the one with ampli- 
tude K2, and we are lead to believe that the star has a planet 
companion with about one Neptune mass. 

In the fourth example the star is being perturbed by a close 
hierarchical binary composed of a star with 1 Mq and another 
body with 0.1 Mq. We have K\ at the current resolution level 
while K2 is 5 times smaller. Again we are lead to believe that the 
star has a planet companion with 1.4 Neptune masses. 

In the next two examples the binaries are composed of a star 
with 1 Mq and a brown-dwarf with 10 Mj (Ex. #5) or a planet 
with 1 Mj (Ex.#6). In these cases both peaks are currently unde- 
tectable (~ cm/s) but could be within reach of fu ture observation 
campaingns like CODEX dPasquini et alj|2008l) . 

Finally, in the last two examples we have a 1 Mq star and 
a binary composed of two Jupiter-sized planets (Ex.#7) and a 
Jupiter-sized planet and an Earth-sized "satellite" (Ex.#8). In 
both cases, the inferred values for M p are too low for planets. 



5. Discussion and conclusion 

In this paper we studied a triple system composed of a star which 
we observe and a binary system. Our goal was to decide if the 
effect of the binary system on the star's center of mass motion 

could be mistaken as a plane t companion to it. 

Previous work done by Schneider & Cabrera (2006) sug- 
gested that the effect of a star-star system could mimic the pres- 
ence of a planet companion to the observed star. However, we 



8 



M.H.M. Morais and A.C.M. Correia: Stellar wobble caused by a binary system 



found that this is not correct, since a more careful analysis of 
the observational data should allow us to distinguish between 
the two effects. Indeed, when there is a close binary system, the 
Fourier spectrum of the star's radial velocity should exhibit not 
one but two peaks located at nearby frequencies. This could also 
indicate the presence of two planets on close orbits, but since 
these systems are most likely unstable we should be able to re- 
ject this hypothesis. 

Nevertheless, we saw that it is sometimes possible to mistake 
(at least temporarily) the effect of a binary system by a planet. 
This happens when the observations are near the resolution 
limit so that we can identify the primary peak but not the sec- 
ondary peak. However, even in this case ou r estimate for the fake 
planet 's mass differs from that obtained bv lSchneider & Ca brera 
(120061) . 

In this article we obtained tools for distinguishing between 
the two effects: we estimate the parameters of the binary system 
that could me mistaken by that planet and check if they are real- 
istic. We recommend that, in future discoveries of planets around 
stars in binary systems, these formulas should be used in order 
to exclude the possibility of a false alarm. 

We also wanted to evaluate the effect of different types of bi- 
nary systems on the star's CM motion. In particular, we were in- 
terested in star-star, star-planet and planet-satellite systems. We 
saw that the effect of a star-star system is only significant when 
this is close to the observed star. For instance, a binary system of 
dwarf stars at 1 AU from each other and at the distance 10 AU 
from the observed star causes a wobble on the latter which would 
be at the current resolution level. We also saw that the effect of 
a star-planet system is currently undetectable although it is pos- 
sibly within reach of future observation campaigns. Finally, we 
saw that the effect of a planet-satellite system is insignificant and 
should never be mistaken by a planet companion to the observed 
star as the implied mass would simply be too low. 

Our conclusions were established for a triple system of bod- 
ies evolving in co-planar and circular orbits. However, it is 
known that these kind of systems may present eccentric orbits 
and different mutual inclinations. The inclusion of these two ef- 
fects may change some of the results presented here. 



Appendix A: 

We re-write the equation of motion of the star M* (Eq[T9ll using 
Cartesian coordinates in the (i, j) frame (see Fig. 1) such that 



r 

u = 
with 



cos 0i + sin9j , 
- sin9i + cos 9j , 



e = e + nt, 

where 9o is an initial phase (FigfJJl. 
Thus we obtain 



X = — 



GM 



1 + -p 2 — Icos^ 
4 F Ml 



+ —p 2 — cos 25 cos 9 + —p 2 — sin 2S sin ( 
4 M 2 M 



Y = — 



GM 



1 + -p 2 — \sin( 
4 F Ml 



(A.l) 
(A.2) 



(A.3) 



(A.4) 



(A.5) 



^ 2 A* ^ 2 ^ -ir 

+ —p — cos 25 sin 9 p — sin 25 cos < 

4 M 2 M 



where p - ai,/a <K 1. 

In order to solve this set of differential equations 
the transformation of coordinates 



X \ _ I sin 6 - cos 0\( X - a* cos ( 
Y ]~ cos0 sin6> Y - a*sin6> I ' 



with 



M (. 3 2 fi 

1 H — p — a 

M + M+X 4 M ] 



we perform 
(A.6) 

(A.7) 



The set of differential equations becomes 

X-2QY-Q 2 X = -l2aj 2 6sin2S , 

Y + 2QX-S1 2 ? = -18w 2 (5cos2S , 

with 

8M\a ) 

This admits the solution 

X - 6x&in2S , 

Y - 6ycos2S , 

where 

- 4(w - Clfdx + 4Q(w - Q.)6 Y - Cl 2 5 x = -\2ur5 
-4(w - Q) 2 5 Y + 4Q(w - Q.)5 X - Q 2 5 Y = -\%u 2 6 

The solution to Eqs. flA.13t and flA.14t is then 

m-f3) 2 +{3 2 + 6/3(1-/3) 



6 X = 126 
6 Y = 126 



(4(l-/3) 2 -/3 2 ) 2 
6(1 -p) 2 + 3/3 2 /2+ 4/3(1 -/3) 



(A{\-/3) 2 -/3 2 ) 



■2\2 



with/? = Q/cj <k 1. 
To zero order in (3, 

6x = 36, 
6 Y = 4.55, 

hence 



a* cos ( 
a* sin t 



sinf 



cos I 



-cos 6 sin 9 



3 6 sin 2S 
4.5 6 cos 25 



(A.8) 
(A.9) 



(A. 10) 



(A. 11) 

(A. 12) 



(A.13) 
(A. 14) 



(A.15) 
(A. 16) 



(A.17) 
(A. 18) 



(A. 19) 



Note that this does no t match the solution obtained in 
ISchneider & Cabrera! (120061) : 



- The 2nd term in Eq.( |A.19t is rotated by an angle 9. 

- The frequency of the 2nd term is 2(w-Q) « 2a> as we assume 
Q. « a). 

- The expressions for Q. (EqIT4Tl and a* (Eq |A.7t take into ac- 
count the mass M+. 
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